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I. INTRODUCTION 

QCD factorization is the fundamental concept to provide theoretical grounds for applying Perturbative QCD to the 
description of hadronic reactions. According to factorization, any scattering amplitude A in QCD can be represented 
as a convolution of perturbative (E) and non-perturbative (T) contributions. In particular, the representation of the 
amplitude A for the elastic Compton scattering 

j(q) + h(p) -> j(q) + h{p) (1) 

off a hadron target is 

A = J2 E r® T r, (2) 

r 

where r denotes the type of intermediate partons in the convolutions: quarks or gluons. Similarly, the amplitudes Ah 
of 2 — > 2 -hadronic reactions 

hi(pi) + h 2 {p2) -> /i'iOI) + K(P2) ( 3 ) 
can be represented in a general form through more involved convolutions: 

A h =Y,T r ®dg rt) ®T r . (4) 

rr' 

where T and T are parton distributions and A^, is the perturbative amplitude, with the purely gluonic contribution 
A^gg Tt ^ dominating in forward kinematics at very high energies. 

There are two kinds of QCD factorization in the literature: Collinear factorization [J 0] and kx- factorization Q. For 
instance, the DIS structure functions f(x, Q 2 ) are respectively represented in those factorizations as follows: 



f( X , Q 2 ) = j d A^rt) [x/ ^ QZ/fj?)^ ( 5 ) 

and 

f(x,Q 2 ) = [jj ^f ipert Hx/f3,Q 2 /kiMP,kl) (6) 

where f(P ert ) stands for the perturbative components of the structure functions both in collinear and A^-factorization 
while 4> and $ denote the involved parton distributions. In the present paper we discuss only general features of these 
perturbative components, so throughout the paper we will keep for them the same generic notation f(P ert ) in both 
collinear and fc^-factorization. In contrast, we are going to discuss details of the involved parton distributions, so, 
in order to avoid any misunderstanding, we have used in Eqs. (|5I6[) different notations for the parton distributions 
in collinear and fcy-factorization, <j> and $, respectively. Besides collinear and hp- factorization, we will introduce 
in Eq. (| 13[) a more general factorization which we call basic factorization. This factorization involves new, totally 
unintcgrated parton distributions which we will denote ^. The parameter /i in Eq. ([5} denotes the factorization scale of 
collinear factorization. Originally, this scale parameter was introduced a kind of border between the perturbative and 
non-perturbative domains of QCD. In addition, it plays the role of the cut-off for the infrared-divergent perturbative 
contributions. Also, it is often associated with the starting point of Q 2 -evolution. Collinear factorization Eq. ([2]) treats 
the intermediate partons as nearly on-shell ones. In particular, this takes place in DGLAP0] and its generalizations 
(see e.g. the overview |5j) to the small-x region. The collinear factorization cannot be used when the perturbative 
contributions to Eq. are calculated using BFKL@, where the external gluons are kept essentially off-shell. In order 
to embrace this case, collinear factorization was replaced by kx -factorization, kr- and collinear factorization, being 
introduced for such different motivations, look unrelated to each other. In Ref. Q we showed that both factorizations 
for the DIS structure functions can be derived from a more general factorization, which we addressed as the Basic 
form of QCD factorization, where the convolutions are totally unintcgrated. 
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In the present paper we continue to investigate the relations between fey- and collinear factorization in more detail 
than was done in Ref . Q , paying more attention to the gauge invariance of the amplitudes involved and considering a 
more complicated dependence of the parton distributions on the transverse momenta. We first consider DIS structure 
functions and then generalize the results obtained for them to the more intricate convolutions describing factorization 
of the amplitudes of 2 — > 2-hadron scattering. Throughout the paper we focus on the high-energy (small- a;) domain. 
Our paper is organized as follows: In Sect. [TT] we consider the derivation of the collinear and fey-factorization for the 
amplitudes of forward Compton scattering off a hadron target and then proceed to the amplitudes of the forward 
hadron scattering in Sect. IIII1 In Sect. HVI we compare the form of collinear factorization obtained in Sects. ITTllIIII with 
the conventional form of this factorization and relate these two forms. In Sect. |V]we discuss the ultraviolet behavior 
of the totally unintcgratcd parton distributions and formulate the requirements for the factorization convolutions to 
be UV-stable; applying these requirements to the standard fits in Sect. I VII we argue against the use the singular 
factors x~ a in the fits. We summarize our results in Sect. IVII1 



II. FACTORIZATION FOR FORWARD COMPTON SCATTERING 

The scattering amplitude A^ v (p, q) of forward Compton scattering off a hadron target is depicted in Fig. [1] 
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FIG. 1. Forward Compton scattering off a hadron target. 

where the blob signifies the presence of perturbative and non-perturbative QCD contributions. We have used the 
standard notations: p stands for the hadron momentum and q (—q 2 = Q 2 ) is the virtual photon momentum. In what 
follows we omit dependence of on the hadron mass and spin. Due to the Optical theorem, the imaginary part 
is proportional to the hadronic tensor W^ v . Using the standard projection operators, A M „ can be expanded 
into a set of invariant amplitudes A r , every DIS structure function f r being expressed through a certain A r : 



f r = -SA r . (7) 
7T 

The next step is to represent A r in a factorized form as an infinite set of i-channel convolutions, each with a certain 
number of the intermediate i-channel partons (quarks and gluons) as shown in Fig. [5] Each of the blobs in Fig. [5] can 
contain both perturbative and non-perturbative contributions. In the present paper we consider only the simplest 
convolutions (a) and (b) in Fig. [2] with two-quark and two-gluon intermediate states respectively, denoting them by 

Ar and A r respectively. In other words, we assume that 



A r ^A r ^+A<fK (8) 

Let us notice that our assumption is in full agreement with the conventional approach to factorization. On the 
other hand, neglecting the contributions with multi-parton intermediate states could rise questions about the gauge 
mvariance of A ( r q) and A { r a) . We discuss the gauge invariance of them in the present Sect. We will focus on discussing 
A r q ^ while the convolution A r 9 ^ with two-gluon intermediate state can be studied quite similarly (see Ref. Q for 
detail). In the analytical form, A r q ^ is 



4 



4 q) (p,q) 



d 4 k kA r (q,k)kT r (p 7 k) 



(2tt) 4 k 2 k 2 

where A r (T r ) denotes the upper (lowest) blob in Fig. [2fa) and k stands for the intermediate quark momentum. 



(9) 
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FIG. 2. Representation of Compton amplitude through convolutions. 

Each of A r and T r can contain both perturbative and non-perturbative contributions, so the factorized form of 
A r in Eq. ([9]) does not correspond to the concept of QCD factorization where the perturbative and non-perturbative 
contributions should be separated. Because of this reason we refer to Eq. (j9|) as the Primordial Convolution. The 
integration in Eq. © runs over the whole phase space and includes both the ultraviolet (UV) and infrared (IR) 
regions. The perturbative contributions in A r are IR-divergent and must be regulated. In the first place, there is 
the IR-divergent power contribution 2pk/k 2 . Such a contribution appears in one of the amplitudes A r with vacuum 
quantum numbers in the t-channel. Throughout the paper we will refer to such an amplitude as the singlet amplitude 
and denote it by As, with the upper blob As and the lower blob Ts, omitting the superscript q. Through the Optical 
theorem As is related to the singlet DIS structure function F\ . We will refer to all other amplitude as non-singlet ones 
and generically denote them Ajvs, with the upper and lower blobs A^s and T^s respectively. We stress that some 
of such non-singlet amplitudes are related to the flavor singlet structure functions (for example, F2,gf, etc). Besides 
the power IR-dependent term, there are IR-divergent logarithmic contributions ~ \n n (2pk/k 2 ) and ln"(Q 2 /fc 2 ). Such 
contributions exist in the singlet and non-singlet amplitudes. It was shown in Ref. [7] that in order to regulate the IR 
divergences in As and Ans without explicitly cutting off the region of small k 2 from the integration region in Eq. (|9|). 
the amplitudes Tg,T^s should obey 

T NS ~(k 2 y, T s ~(fc 2 ) 1+7 (10) 
at small k 2 , with 7 > 0. It is convenient to rewrite Eq. (|9|) in terms of the standard Sudakov variables 



k = -a{q + xp) + (3(p — x'q) + k±, 



where x = —q 2 /w, x' — p 2 /w and w = 2pq. Doing so, we arrive at 



kA r (q, k)kT r (p, k) 



(wa/3 



(11) 



(12) 



where r — S, NS. Now let us consider A^ (p, q) in the Born approximation. The Compton amplitude off the quark 
(the upper blob) in the Born approximation is depicted in Fig. [3] 

The quark amplitudes A$ NS and the gluon amplitudes A^ in the Born approximation were investigated in detail in 
Ref. Q, so we omit this in the present paper. Let us stress that the lower blobs Ts and TVs' in t ne Born approximation 
are altogether non-perturbative. As is known, the Born Compton amplitude is gauge-invariant when the quark is 
on-shell. In Appendix A we demonstrate that gauge invariance for those amplitude is restored in the high-energy 
region where w Q 2 ,k 2 . 
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FIG. 3. The Born contribution to Compton amplitude. 
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FIG. 4. Examples of the graphs contributing to Compton amplitude beyond the Born approximation. 

Going beyond the Born approximation means adding extra gluon and quark propagators in all possible ways to the 
Born graphs. Examples of this are shown in Fig. |4] 
The resulting graphs can be classified into two groups: 

(i) The graphs where extra propagators do not affect the lower blob (for instance the graphs (a) and (b) in Fig.0}. 

(ii) The graphs where the extra propagators involve the lower blob (graph (c) in Fig. [4]). 

Obviously, accounting for the whole set of graphs from group (i) leads to the convolution depicted in Fig. [U[a). 
whereas graphs from the group (ii) form convolutions with more complicated intermediate states. Such graphs cannot 
be simply neglected because this affects the gauge invariance. However, we demonstrate in Appendix A that gauge 
invariance is restored at small x even when the group (ii) is neglected. Let us stress that we have arrived at the 
convolution where the upper blob A r is altogether perturbative and at the same time the lower blob T r is totally 
non-perturbative. 

Applying the Optical theorem to Eq. (|I2p . simplifying the spinor structure and adding the similar gluon contribution 
(see Ref. [7[ for detail) leads to the convolution for the DIS structure functions: 



fs, NS (x,Q 2 ) = I d<x% f 'I (s/A Q 2 /k 2 )* s ,Ns(wa, fc 2 ) (13) 

J p (wap + Kj_) ' 

where f^g"^ {fs'^) is the perturbative component of /ns (fs) and ^ns — ST/vs, = QTs- In order to regulate 
the IR-divergences in Eq. (fl3|) . these parton distributions should decrease with k 2 at small fc 2 : 



* NS ~(k 2 )\ * s ~(k 2 ) 1+1 - (14) 

Eq. (|13[) represents the structure functions in factorized form, with perturbative and non-perturbative contributions 
being separated, so this representation authentically corresponds to the concept of QCD factorization. On the other 
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hand, this factorization involves three integrations, which distinguishes it from both kx -factorization (which involves 
integration over /3 and k±) and the collinear factorization (involving integration over /3only). InRef. we introduced 
the term " Basic form of factorization" for such totally unintegrated convolutions and will use it throughout the present 
paper. The basic factorization involves the totally unintegrated parton distributions. In this respect we are close to 
Rcf. [9j but in contrast to that paper, we focus on small- a; kinematics and go beyond the Born approximation. 
The basic factorization can be approximately reduced to the kx -factorization if the restriction 



wa/3 < k\ (15) 

is accepted. In this case one can perform the integration over a without dealing with f(P ert ) and arrive at the kx 
-factorization for fg, /ns'- 



fs,Ns(x,Q 2 ) « fj£ /S(^^ 2 A1) ^sM^kl), (16) 

$s and §ns being the parton distributions for the kx -factorization. They are related to the totally unintegrated 
parton distributions ^g, ^jvs 1 

$jvs = / da^ NS (wa, k 2 ), $ s = / da^ s (wa, k 2 ). (17) 

I (0) I I (°) I 

The lowest limit of integrations in Eq. (|17[) is fixed from the following considerations: the invariant energy S2 = (p — k) 2 
of $ns an d $s must be positive, so 

s 2 = (p- k) 2 &wa-k 2 ± > 4 0) > 0. (18) 
In order to regulate the IR divergences in the fc^-convolutions, these parton distributions should behave (cf. (fTU|) ) 

as 

* JVS ~(*i) 7 , $ s ^(fci) 1+7 (19) 

at small k 2 ^. Obviously, the fcT-factorization in Eq. (|16p authentically coincides with the conventional form of fc-r- 
factorization in Eq. ©. 

In its turn, fcx -factorization can also be approximately reduced to collinear factorization. In order to do it, we 
assume the peaked dependence of the parton distributions on k"]_ as shown in Fig. [5] 

We stress that this hypothesis can be checked by analysis of experimental data. The peaked dependence $s\ns on 
k\ makes it possible to perform the integrations over k\ in Eq. (|16[) without involving the perturbative components 

fgjfg , arriving at the collinear factorization 



fs,Ns(x,Q 2 ) « ]T f ^ftTs^/^Q 2 /^) <f>S,NS(P^l), (20) 

k Jx P 

with the singlet and non-singlet parton distributions 4>ns{P, H 2 ) and ^s(/3,/i^) defined at the factorization scale /j.^. 
They are related to the kx -parton distributions in the following way: 



'JNS — 



dk\^ NS {wa,k]_), 4> 



S = 



dk^siwa^k 2 ^). 



(21) 



The integration regions Dk in Eq. (12 1 1) are located around the positions of the maxima of the parton distributions at 
k\ = fj, 2 , as shown in Fig. [SJ Obviously, Eq. (1201) differs in many respects from the conventional expression in Eq. ([5]) 
for the DIS structure functions in collinear factorization. We will show how to bring Eq. (|2"U|) to the conventional 
form in Sect. IV. 
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FIG. 5. Assumed dependence of the parton distributions on k\ 
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FIG. 6. Forward 2 — > 2 scattering of hadrons. 

III. FORWARD SCATTERING OF HADRONS 

Let us consider the 2 — > 2-scattering of hadrons (see Eq. ([3])) in forward kinematics where s = (pi + p 2 ) 2 ^ t = 
{p'l — Pi) 2 - This process is depicted in Fig. [HI 

In order to be able to apply Perturbative QCD to this reaction, should be regarded as consisting of the following 
two sub-processes as depicted in Fig. [7J 

(i) Each of the colliding hadrons emits one or several partons (the upper and lowest blobs in Fig. [7]). 
(ii) High-energy interaction of the emitted partons (the middle blob in Fig.[7J. 

In the present paper we will consider convolutions with two intermediate partons in the i-channel (graphs (a) and 
(b)) only. This approximation allows us to write the scattering amplitude Ah of the process of Eq. ([3]) in the form of 
the following convolutions (cf. Eq. @): 

A h w A "- = E T (r ' } (k 2 , k' 2 , P2 , P ' 2 ) ® H^ r ' r \k 2 ,k' 2 , ki, k[) T^(p 1 ,p' 1 ,h, k[), (22) 

rr' rr' 

where the superscripts r, r' refer to the kind of intermediate two-parton states: two-quark or two-gluon states. and 
T^ r ^ stand for the lowest and upper blobs respectively while H^ r r ' denotes the middle blob in Fig. [Jj k' 12 = q + ki_ 2l 
with q = p'i — pi ■ 

Similarly to the Compton scattering amplitudes considered in Sect. II, the amplitudes Ah for hadron scattering can 
have either vacuum or non- vacuum quantum numbers in the i-channel. In the present paper we focus on the singlet 
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FIG. 7. Representation of the amplitude Ah through convolutions of sub-processes. 



amplitudes with the vacuum quantum numbers. Then, in what follows we consider in more detail the amplitude A, 
which corresponds to graph (a) in Fig. [7j 



A aa = J ^^ U^^^Lfe^^Uti,*!). (23) 

where the subscripts fi, (jf, A, A' denote the gluon polarizations. In order to account for the leading contributions, we 
parameterize T? 9 ?, and TjfJ, as follows: 



j,(g) _ 2p 2M P2 M ' j, T (g) _ ^PlXPlX' j, „ _ 2p2 M P2 M ' „(gg) ^V\\V\\' 

fill 1 s ' AA' s ' s fifi'W' s ' y ) 

T, T, H being scalar functions. We have dropped the superscripts here. Using these notations and the Sudakov 
parametrization of Eq. (|lip for ki.2, we can rewrite Eq. (1231) as follows: 



A gg = J da^d^dkl^fis^klk^^M^Tisauklki 2 ), (25) 

with k 2 = —sa r (3 r — kf ± . Let us explain the notations we have introduced in Eq. (|25|) . We have denoted by s' the 
invariant energy for the sub-process of gluon 2 — > 2 -scattering: s' = (ki — k-i) 2 = (k[ — k' 2 ) 2 ~ so^i- Throughout 
the paper we consider the kinematics where 

s 1 ->klklk' 2 ,k 12 . (26) 

Let us denote by k c the momentum with the largest virtuality: k 2 ~max[k 2 , k 2 , k' 2 , k' 2 }. Then we denote by k 2 
any of k 2 ,k' 2 , providing b ^ c. Similarly, k 2 is any of kf,k' 2 , providing a ^ c. Finally, k 2 —mm[kf,k2,k' 2 ,k' 2 }. 
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The amplitudes T, T can include both perturbative and non-perturbative contributions, as can amplitude M. The 
amplitude M replaces the amplitude H when the momenta k ai kh,k are used instead of hi, k%. The integration over 
fci ; 2 in Eq. ()25[) runs over the whole phase space, so the integration region includes the IR-singularities where any 
of k 2 , k' 2 or fcf, k 12 can be equal to zero. Obviously, fcf and k' 2 2 (and k\ and k' 2 ) can be small at the same time 
only when q — 0, i.e. when t — and therefore the forward kinematics becomes collinear. This kinematics is the most 
IR-singular. At t = Eq. (|25|) looks as follows: 



A gg = J da h2 df3 h2 dkl 2± T(sf3 2 ,k 2 2 )^ (j^j M^T{sa x ,kl) (27) 

where k 2 = mm[k 2 ,k 2 ]. 

The perturbative infrared-sensitive contributions in the amplitude M are logarithms: M = M (\n(s' /kf), ln(s'/ /cf )) . 
On the other hand, the hadron amplitude A gg must be free of IR problems by definition. In order to keep Eq. (f2"7| 
IR-stable, Amplitudes T and T should compensate the IR-divergent terms (cf. Eq. ([TD])): 



T~{kl) 1+ \ T ~ {k 2 ) 1+1 , (28) 
with 7 > 0. When t ^ 0, the conditions for the IR stability are weaker than in Eq. 



T~(*£r, T~(k 2 b y. (29) 

With the integrations in Eqs. (|25|27p IR-stable and the total energy of the gluon scattering high, the amplitude M 
becomes completely perturbative. When M is calculated in the Born approximation, the amplitude Ah is represented 
in Fig. [S] Let us notice that the Born approximation means that the blobs in Fig. [5] include non-perturbative 
contributions only. 



FIG. 8. Convolution for Ah in the Born approximation. 



Considering A gg beyond the Born approximation corresponds to adding the quark and gluon propagators to the 
Born graph. In doing so, we account for those graphs which do not involve the upper and lower blobs, like graph (a) 
in Fig. IHl but avoid including the graphs involving the blobs like graph (b) . 

Such a procedure keeps the blobs totally non-perturbative. The gauge invariance of the set of retained graphs 
is proved in Appendix B. As a result, we arrive back at graph (a) depicted in Fig. [7] where, however, the middle 
blob is perturbative while the upper and lowest blobs are altogether non-perturbative. Therefore, the convolutions 
in Eqs. (I23I25P complemented by the IR-regulators of Eqs. (|2"51 |2U)) represent the hadron scattering amplitude Ah in 
terms of the convolutions of Eqs. (|25l27p . with the convolutions being in the perfect agreement with the concept of 
the QCD factorization. Applying the same analysis to the amplitudes A qq , A qg , A gq , we arrive to the representation 
of the hadron amplitude Ah in terms of the QCD factorization. Obviously, this factorization differs from the fcx 1 - and 
collinear factorization, so we again refer to it as the Basic Factorization. 
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(b) 



FIG. 9. Convolutions for Ah beyond the Born approximation. 



A. Transition to the kr -factorization 



The transition to kr -factorization means performing the integrations over a\ and f3 2 in Eqs. (I25I27P without dealing 
with the amplitude (s'/fc 2 )M. This is possible only if these integrations are restricted as follows: 



2±- 



(30) 



The meaning of these restrictions is that the virtualities of the intermediate partons are space-like. This agrees with 
the known fact that the virtualities of the ladder partons in the Leading Logarithmic Approximation are space-like. 

Restriction (|30|) allows us to integrate over fa and a± in Eqs. (|25l27p . arriving at the following expression for Ah(s, t) 
in fcx -factorization: 



A h (s,t) = I da 2 d/3 1 dfc^dfc 2 \*( S a 2 ,fc 2 \,fc^)^A( pert ^*( S /3 1 ,fc^,fc;i), 
with the perturbative amplitude A^ pert > and non-perturbative parton distributions "J, VP defined as follows: 



(31) 



A ( P ert) = ( s 7^)M(ln( S 7fc 1 2 ^),ln( S 7fci),ln( S 7fci 2 ± ),ln( S 7^ 2 )) , 



(32) 



*= / dfaT(sfa,kl ± ,k' 2 \), * =/ daiT(a ai , fcg.) 



(33) 



with Si and s 2 being the minimal invariant energy of the upper and lowest blobs, respectively (cf. Eq. (|18[1 ). 
Putting t = 0, then taking the imaginary part of Eq. (I31[) and using the Optical theorem, we arrive at the expression 
for the total cross-section a to t of the hadron scattering in kx -factorization: 



„,„ = / ,/,, ,1/3, f J-l fl-L $(aa 2 , fc 2 2 ± ) tr£p> (*', fc 2 ± , fcfj fc 2 J, 



'hi. "^J. 



(34) 



with the new parton distributions $(s/?i, A;^) and $(sa 2 , fc|j_) defined quite similarly to the parton distributions in 
Eq. (TTTD: 
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$(s/3 1; fc 2 J = fc? ± ), $H, fc 2 \) = 9*(«a 2 , k 2 ± ). (35) 

In Eq. we have used the notation af^^ for the total cross-section of the perturbative sub-process describing the 
scattering of two partons. fcr-factorization was introduced in Ref. Q to replace the collinear factorization in the case 

when the BFKL PomeronQ is used for the description of <jf^ . Collinear factorization can be used when alternative 

methods are exploited to describe a^ rt \ 



B. Transition to collinear factorization 

The transition from fc^-factorization in Eq. (|34p to collinear factorization can be done, assuming a peaked k±- 
dependence of the parton distributions involved, $, $. This dependence is depicted in Fig. [SJ Once we accept this 
assumption, we can neglect the k±- dependence of c^tot^ an< ^ integrate only the parton distributions <J>, $ over fci.2j_- 
As a result, we arrive at 

vtot =J2 da ^ d Pi 9k(sa 2 , Hi) cr^ t ert) (s',/4) ¥> fe (s/3i,/4), (36) 

k,l J 

where 



(p k {sp 1 ,fil)= dkl^isPx,^), (p k (sa 2 ,(j.l)= dk?; ± <S>(sa 2 ,k 2 1± ) (37) 

JD k JD k 

where D k is the region around the maximum of the parton distributions at k 2 2J _ = 

The transitions from the basic form of QCD factorization to kr -factorization and to collinear factorization for the 
non-singlet component of Ah can be done quite similarly. Obviously, Eq. (|36l) differs from the conventional expression 

a to t = J da 2 d/3i(j} k (sa 2 ,^ 2 ) cr^f*' (s', ^ 2 ) 0(s/?i,^ 2 ), (38) 
where the parton distributions </>, <fi include both perturbative and non-perturbative contributions. 



IV. COMPARISON OF EQS. (|20|34|) TO THE CONVENTIONAL FORM OF COLLINEAR 

FACTORIZATION 

The reduction of the kx -factorization to collinear factorization led us to Eq. ([2Hf for the DIS structure functions 
and Eq. (|36[) for the total cross-section. Obviously, these expressions differ a lot from the expressions written in 
the conventional form Eq. (15I38[) of collinear factorization. Let us show how Eqs. (j2QI36[) can be brought to the 
conventional form. First of all, notice that the collinear factorization in Eqs. (120136)) differs from the conventional 
form in the following aspects: 

(A) There is only one factorization scale /i in Eqs. (I5I38[) and and this scale does not bear any physical meaning, 
so /i can be chosen arbitrary. On the contrary, Eqs. (|2QI36p admit the possibility of having several intrinsic scales /i r , 
each corresponds to the maximum k\ = ji 2 , of the unintcgrated parton distributions. 

(B) The parton distributions ip, <p in Eqs. (|20I36[) include the non-perturbative distributions only whereas the 
conventional form operates with the distributions (f>, 4> where there are both perturbative and non-perturbative 
contributions. 

Despite such a considerable difference between these two approaches, it is easy to show that in fact Eqs. (120136)) 
can be brought to the form of Eqs. (|5l38p . This can be done identically for Eq. (|2T)]) and Eq. (|3"o| . so in what follows 
we will focus on considering Eq. (|20|) . Obviously, Eq. (|5|) can be written as (cf. Eq. ((2J) 



f(x,Q 2 ) = E(Q 2 /v 2 )®j>(p 2 ) 



(39) 
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where E(Q 2 //j J 2 ) denotes any of the appropriate evolution operator for performing the evolution from u 2 to Q 2 . We 
have dropped its x-dependence as unessential. In particular, it can be the DGLAP -evolution operator[4|. Using such 
notation, Eq. (|20|) can be rewritten as follows: 



f(x, Q 2 ) = ]T E(Q 2 /^ 2 ) g fafal) = E(Q 2 /[i 2 ) ® E(^/f4) ® Mti) = E(Q 2 /fi 2 ) ® <j>(p 2 ), (40) 
fe fe 

with 

<j>(» 2 )=J2 E (vV»l)®Mf4)- (41) 
fc 

Eq. (j4"0l makes sense if [i> [ik- This perfectly agrees with the actual situation: the scale /x in the standard approach 
is usually chosen ~ lGeV whereas the intrinsic scales fik are supposedly non-pcrturbative and they are therefore 
expected to be not far from Aqcd- Eq. (|40j) demonstrates that the conventional approach and our approaches can 
be converted each into other when the factorization scale \x is chosen at the GeV range, so they are equivalent and 
indistinguishable when \x is kept within the GeV range or higher. In this case the parton distributions should explicitly 
include perturbative terms in addition to the non-perturbative ones. They become distinguishable at smaller values of 
the factorization scale, close to Aqcd (where the perturbative component vanishes), providing there are several (more 
than one) intrinsic scales /i^, otherwise they are again indistinguishable. However, the difference between them can 
be found in the framework of the hr -factorization where the existence of the maximums in dependence $ = 
can be checked with analysis of experimental data. 

V. ULTRAVIOLET BEHAVIOR OF THE PARTON DISTRIBUTIONS T s , T NS 

The basic factorization convolutions for the Compton amplitudes (i.e. Eqs. (|9ll2p complemented by Eqs. (fTU)) ) and 
for the hadron amplitudes Ah (i.e. Eqs. (|25I27|) complemented by Eqs. (|28|) ) are IR-stable. In addition, they should 
be UV-stable. The UV-stability of the Compton amplitudes was discussed in detail in Ref. Q- It was proved that 
the convolutions are UV-stable when the parton distributions Ts (a) , T/vs (ct) at large |a| behave as follows: 

T NS ~ la]- 1 -**, T S ~ \a\-\ (42) 

with h > 0. It is easy to show that the singlet parton distributions T(a>i), T{p2) in the convolutions of Eqs. (|25|27p 
should exhibit absolutely the same UV behavior: 

T( ai )~\ ai \- h , f~^ h . (43) 

VI. RESTRICTIONS ON THE PARTON DISTRIBUTIONS 

We obtain below the restriction on the parton distributions in the collinear and fcr-factorizations that follow from 
the integrability of the basic convolutions. We exploit here the obvious mathematical requirement of stability of the 
convolutions in both the IR and UV -regions. These restrictions can also be used as suggestions when PDF fits are 
constructed. 

A. Parton distributions in kr -factorization 

The fits for the parton distributions in fcT-factorization are commonly fixed from phenomenological considerations 



(see Ref. [ll|]). Below we give some restrictions on the fits following from theoretical grounds. The parton distributions 
$ in fc^-factorization are defined in Eq. (jTTJ) (see also Eq. (|3"5|l ) as the integrals of the totally unintegrated parton 
densities. Combining this equation with the restriction of Eq. (fTUj) (see also Eq. (|28p) on their IR-behavior leads us 
to the most general form of $ns, $s'- 

$ns = (fci ) 7 x^(/?. k\) + (k 2 ± y- h /} h xP s (P> fc i). ( 44 ) 
$5 = (fci) 1+ ^ 1} (/3, fci) + {kl?^- h r 1+h x { s\Pi hi), 
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with the requirement 



7-ft>0 



(45) 



We remind that the positive parameters 7 and h were introduced in Eqs. (|10I42I) for the Compton amplitudes to 
guarantee their IR and UV stability, respectively. When the Compton amplitudes are considered in basic factorization, 
7 and h are independent. However, they proved to be related by Eq. (|45|) when basic factorization has been reduced 
to kr -factorization: Eq. (|4"5|) ensures the IR-stability of the structure functions in fc*r-factorization. Eq. (|4"4l shows 
that the parton densities in fc^-factorization consist of two terms, each involves different power factors of kj_ and the 
last fc^-factor involves the longitudinal variable j3. In addition, they involve arbitrary functions XsnS' ^ u ^s 'ns 
const at k\ — > 0. $s,NS originate from Ts.ns, so the requirements of the UV stability in Eqs. (|42I43I) mean that the 
fits for $ns should not contain the factors f3~ a , with a > 0, whereas the requirement on fits for <I>s is less severe: they 
can contain such factors, providing a < 1. However, sometimes the same parton distributions contribute to the singlet 
and non-singlet constructions. For example, this is true for the flavor singlet components of the structure functions 
Fi (addressed as the singlet in the present paper) and F 2 (addressed as the non-singlet): they both involve the same 

parton distributions. In such cases, the singlet fits should obey the more severe requirement for the non-singlets. 

(12) 

In order to make possible the transitions to the collinear factorization, the fcx-dependence of Xs jvs snou ld be of a 
peaked (for example, Gaussian) form. 



The conventional fits[12j for the parton distributions in collinear factorization are also constructed from purely 
phenomenological considerations. Below we impose some restrictions on these fits arising from integrability the 
factorization convolutions. The standard fits for the initial parton densities <5q, Sg in collinear factorization are known 
to include a normalization N, the singular factors x~ a , with a > 0, and regular terms. For example, 



where the parameters N, a, b, c, d are positive. They are specified from analysis of experimental data. Although such 
expressions do not look explicitly like the ones obtained with the perturbative methods, nonetheless there are two 
options to study: 

(i) The singular factors x~ a appear as a result of total resummations of the perturbative contributions. In this case 
we identify the fits with the parton distributions <fi containing both perturbative and non-perturbative contributions. 
The factors x~ a in such distributions mimic the resummation of either the leading logarithms (see [l(|) or sub- leading 
logarithms (see Q). In the both cases the singular factors can be removed from the fits when the resummation is 
accounted for because the resummation of logarithms of x can be absorbed by the coefficient functions. 

(ii) The whole fits Eq. (|46|) . including the singular factors x~ a , have the non-perturbative origin. In this case we 
identify the fits in Eq. (|4"6l with the non-perturbative distributions tp defined in Eqs. (I21I37[) and apply to them 
the restrictions of Eqs. (|42I43[) . These restrictions exclude the use of the singular factors in the expressions for the 
non-singlet structure functions F2,F^ S \g%, etc and also suppress the singular factors with a > 1 in the expressions 
for the singlet F\ . However, the parton distributions used for F\ and F2 are identical, therefore the suppression of the 
singular factors with a > can be applied to all structure functions, including the singlet F\. The singular factors 
x~ a in the DGLAP fits for initial parton densities should be removed from the fits because they contradict to the 
integrability of the basic convolutions of the Compton amplitudes. Let us remark that the removal of the factors x~ a 
from the fits and replacing such factors by the total resummation of the logarithms can reduce the fits in Eq. (|4"o| 
down to much simpler expressions like N(l — x) b or, when the starting values of x are small, down to constants. 



We have shown that the scattering amplitudes of the processes Eq. ([T]) and Eq. Q can be represented as basic 
convolutions where the perturbative and non-perturbative contributions are separated (i.e. located in different blobs). 
The basic convolutions correspond to concept of QCD factorization but they are more general than collinear and 
fey-factorization. The perturbative components of the basic convolutions are off-shell and therefore they are not 
gauge- invariant. However, their gauge invariance is restored at high energy (small x). When the virtualities of the 



B. Parton distributions in collinear factorization 



5q = Nx- a {l-x) b (l + cx d ), 



(46) 



VII. SUMMARY 
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intermediate partons are space-like, the basic factorization can be approximately reduced to fcr-factorization. In 
its turn, ^-factorization can approximately be reduced to collinear factorization if the peaked dependence on k± 
(see Fig. 5) is assumed for the unintegrated parton distributions and the parton distributions are altogether non- 
perturbative. The sharper the peaks in Fig. 5 are, the higher is the accuracy of the reduction. We stress that this 
hypothesis can be checked by analysis of experimental data. In contrast to the conventional scenario of collinear 
factorization, we have arrived at the more involved form of collinear factorization, with one of several (intrinsic) scales 
fi r corresponding to the maxima of the unintegrated parton distributions. Using perturbative evolution to increase the 
scale, we have brought this complicated many-scale picture to the conventional form, where only one arbitrary scale 
is used explicitly and the parton distributions include both perturbative and non-perturbative contributions. The 
integrations in the basic convolutions run over the whole phase space but they must yield finite results. Exploiting this 
obvious requirement allows us to obtain restrictions on the parton distributions in both collinear and fc^-factorization. 
We have shown that the parton distributions in /cy-factorization should consist of two terms, with different powers 
of fcj_. Then, we exclude the use of factors x~ a from the fits for the parton distributions both in the collinear and 
/cT-factorization. Finally, let us remark that the results obtained in the present paper can easily be extended to 
inelastic hadron scattering at high energy. 
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Appendix A: Gauge invariance of the Compton scattering amplitude 
1. Notation for Compton scattering 

Let us consider the Compton scattering 

l{q) + q(p) ->■ l(q') + «(?')• ( A1 ) 
We denote by A^ v the scattering amplitude of the process (|A1[) . Gauge invariance states that 

q^A^ = q v A» v = 0. (A2) 

2. Compton scattering amplitude in the Born approximation 

In the Born approximation the scattering amplitude A^ is represented by the graphs depicted in Fig. [3] 

A = <p'\7u(p + q + m)j t _ l \p > + <p'h^(p' -g + m)j u \p> 
(p + q) 2 — m? (p 1 — q) 2 — m 2 

3. Gauge invariance of the Compton amplitude in arbitrary kinematics 

It follows from Eq. (|A"3j) that 

q A = <p'\lvq(p-m)\p > _ <p'\(p' -m)q^\p> 
^ (p + q) 2 — m 2 (p' — q) 2 — m 2 

, q 2 + 2pq q 2 - 2p'q 

< p hu\p > 



(p + q) 2 — m 2 (p' — q) 2 — m 2 
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Eq. (|A4[) can be re-written as follows: 



qp,A^ v 



<p'hvq(p-m)\p> <p'\{p' -m)<tf v \p> 



ip + q) 2 — m2 
< v'YiAv > 



p z — m 



ip' ~ l) 2 ~ m2 

/2 2 



(A5) 



(p + f/) 2 — to 2 {p' — g) 2 — m 2 



Obviously, Eq. (|A5j) yields a zero result only when the incoming and outgoing quarks are on-shell. In contrast, the 
photons can be either on-shell or off-shell. 

Let us focus on considering forward Compton scattering. 



4. Forward Compton amplitude in the Born approximation 



In this case 



A,, 



<p\lu{P + q + m)-i^\p> <p\l^{p- q + m)~f v \p> 



(p + qf 



(p — q) 2 — m 2 



(A6) 



Obviously, the forward Compton amplitude is gauge-invariant when the quark is on-shell. However, now we allow it 
to be off-shell. In order to study gauge invariance for an off-shell quark, we rewrite A^ v as follows: 



where 



.4 



- 4(1) + A {2) + A {3) 



(A7) 



A (2) 



A (3) 



<p\lvPlp\p > , < p\lpPlAP > 
[p + q) 2 — m? (p — q) 2 — m 2 

< p\ivqiAp > _ < p\i^qiAp > 
(p + q) 2 ~ m2 {p ~ q) 2 ~ m2 

" <P\lvlp\P> <P\lplv\P > 

ip + q) 2 — m2 (p — q) 2 — m2 

Let us multiply Eq. (|A7j) by q v and represent the result as follows: 



(A8) 



with 



(A9) 



G (i) = A (i) _ 



< p\qpin\p > , < p\lppq\p > 



(p + q) 2 



(p - qf 



(A10) 



^q v A$=q 2 <p\ 1 ,\p> 



(p + q) 2 



(P - q) 2 



(All) 



G^ = q v A$ = m 



<p\tin\p> <p\i^q\p> 



(p + q) 2 — m 2 (p — q) 2 — m 2 



Now let us simplify the Dirac structures in (|A10[) : 



l^qp = -it^XpaqxPplzla + q^P - Ppq + pqip 

qpip = i£n\ P aq\Ppj5ia + qpP - Ppq + pqip 



(A12) 



(A13) 
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This makes it possible to represent G^ as a sum of symmetrical and antisymmetrical contributions: 



(A14) 



with 



G]} A ' = ie fJ ,\ P crq\Pp 



<P\lbJa\p > 



1 — X — Z 1 + X + Z 



(A15) 



where x = —q 2 /w, z = m 2 /w and w = 2pq. 
Similarly, 



G (is) = \% < p\p\p > -pm < pW\p > +pq < p\jv\p >] 



1 — x — Z 1 + X + z 



(A16) 



In order to retain the CP-invariance in (IA15[) . the term < pl^jalp > must yield the spin contribution ~ S a . This 
converts the antisymmetrical factor in (|A15|) into 



^(iXptrQXPpSa- 

In order to simplify (|A16|) we use the approximation 



(A17) 



<php\p >~Pp, 



(A18) 



which allows us to re- write (|A16j) as 



g( p S) ~ [vnP 2 - p^pq + p^pq] 



<P\P > 



1-x 



1 + X + z 



(A19) 



Obviously, the expressions on the r.h.s. of Eqs. (|A15[[A"T9]l are not equal to zero, which means that gauge invariance 
for Ap 1 } is broken. The same is true for A$ and Aj?} when the quark is off-shell. 



5. Restoration of the gauge invariance for forward Compton scattering at small x 

In the high-energy limit where w 3> \q \, \ (p 2 — m 2 )\, i.e. where both 

x < 1 



(A20) 



and 



(is) 



(A21) 



z«l, 

the expression in the square brackets in Eq. (|A19|) becomes small and therefore in this kinematics G : 
0. Then, in this kinematic region the spin S p becomes mostly longitudinal, so the expression in (|A17[) is nearly zero 

0. Eq. (IA11|) demonstrates that Gp in region (IA20[) is proportional to x, so G { p ] » 0. Finally, 



max[i, z\ 



and therefore G 



(1A) 



let us notice that G^ 3 ' ~ m, so it becomes small in the small- a; region (|A20[) . Therefore, the gauge invariance of the 
forward Compton amplitude, although broken when the off-shell quark is off-shell, is restored in the limit of small x 
and z for the unpolarized Compton scattering amplitudes and for those spin-dependent amplitudes, in which the spin 
is longitudinal. 
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6. Accounting for radiative corrections 

The most important radiative corrections in the small- a; region are logarithmic ones. The arguments of the loga- 
rithms can be chosen as w/p 2 and Q 2 jp 2 . Accounting for the radiative corrections with logarithmic accuracy does not 
destroy the Born structure of the forward Compton amplitude. Then, it is obvious that ln(«j/p 2 ) and \a(Q 2 /p 2 ) can be 
large only if p 2 <C Q 2 , w. So, the main contributions in the small- x region come from quark virtualities obeying (|A21I) . 
Therefore the total resummation of such contributions does not break the gauge invariance of the Compton amplitude 
at small x. This gives us the right to consider the factorization convolutions with the two-parton intermediate states 
only and to neglect all other intermediate states with greater numbers of partons without breaking gauge invariance 
at small x. 

Appendix B: Gauge invariance for quark-quark scattering in collinear kinematics 

Let us consider the quark-quark scattering 

q(pi) + ?(pa) -> q(p'i) + q(p' 2 ) (Bi) 

in collinear kinematics where pi w p' ly P2 ~ p' 2 - In what follows we will focus on the colorless, in the i-channel, part 
of the scattering amplitude of this process. In other words, we will consider the amplitude with vacuum quantum 
numbers in the i-channel. In the lowest-order approximation the colorless part of the scattering amplitude Aq is 
depicted in Fig. [10] 

P2 P2 P2 P2 





FIG. 10. Colorless contribution to the quark scattering in the lowest order. 
We write it as follows: 



A = A S + A U , (B2) 



where 

, < Pihuipi ~k + m)7 M k?i >< P2\jv'(p2 + k + m)7^|p 2 > d^^k) d v - v {k) 



((pi — k) 2 — m 2 ) ((p2 + k) 2 — rn 2 ) k 2 k 2 



(B3) 



and 



A ^ <P}HApl ~ k + m^^px >< p 2 \l^'{p2 + k + m)~i v ,\p 2 > d M >(fc) d v , v (k) . 
((pi — k) 2 — to 2 ) {{p2 — k) 2 — to 2 ) k 2 k 2 

As our object is to study the gauge invariance of Aq, we have omitted in Eqs. (IB3IB4|) all factors unessential for that. 
As is well-known, the leading contributions to A s and A u come from the terms ~ Pi,P2 in the numerators of (|B31 
IB4[) , with k and to neglected. The remaining structures are obviously symmetrical under the replacements p ^ v 
and p' ^ v' . Therefore we can write them as 
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(k 2 - 2pik + (pf - to 2 )) (k 2 + 2p 2 k + (p 



»)) k 2 



(B5) 



and 



S^ v {pi)S^ v '(p2) 



dfi'^k) d v > v {k) 



((k 2 - 2 Pl k + ( P 2 - m 2 )) (fc 2 - 2p 2 fc + {p 2 



TO 2 )) fc 2 



A: 2 



(B6) 



with 



S^{pi) =< PihuPi^^Pi >, S^v'ipi) =< P2hn'P27v'\P2 > • (B7) 

Gauge invariance means that the replacement of d u / u (k) by k v (f> v i (fc) , <f> being an arbitrary function, yields a zero 
result. Multiplying S flv (pi) by k v and substituting the result in Eqs. (|B5lB6j) proves the breaking of gauge invariance. 
However, in the high-energy limit where 



|2pifc|, \2p 2 k\ » |fc 2 |,|(p 2 -m 2 )| (B8) 

gauge invariance is restored. This is the appropriate region to produce the leading logarithmic contribution. This 
means that accounting for these contributions does not destroy the gauge invariance of the amplitude Aq. 
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